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Abstract 

A type of adaptive finite element method for the eigenvalue problems is 
proposed based on the multilevel correction scheme. In this method, adaptive 
finite element method to solve eigenvalue problems involves solving associated 
boundary value problems on the adaptive partitions and small scale eigenvalue 
problems on the coarsest partitions. Hence the efficiency of solving eigenvalue 
problems can be improved to be similar to the adaptive finite element method 
for the associated boundary value problems. The convergence and optimal 
complexity is theoretically verified and numerically demonstrated. 
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1 Introduction 

The finite element method is one of the widely used discretization schemes for solving 
eigenvalue problems. The adaptive finite element method (AFEM) is a meaningful 
approach which can generate a sequence of optimal triangulations by refining those 
elements where the errors, as the local error estimators indicate, are relatively large. 
The AFEM is really an effective way to make efficient use of given computational 
resources. Since Babuska and Rheinboldt [5], the AFEM has been an active topic, 
many researchers are attracted to study the AFEM (see, e.g., [2|. f5| f6| fT|. [TOl [T8|. 135] . f36 |. 
H2J HH] and the references cited therein) in the last 30 years. So far, the convergence 
and optimality of the AFEM for boundary value problems has been obtained and 
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understood well (see, e.g., [2H E21 IM1 ESI S31 El EH SOI SH HH H31 HH E] and the 

references cited therein). 

Besides for the boundary value problems, the AFEM is also a very useful and 
efficient way for solving eigenvalue problems (see, e.g., [SJ dSJ E31 EH ESI ESI ED EH])- 
The AFEM for eigenvalue problems has been analyzed in some papers (see, e.g., 
[20| EH EH] and the reference cited therein). Especially, [201 EZ] give an elaborate 
analysis of the convergence and optimality for the adaptive finite element eigenvalue 
computation. In |24j . authors also give the analysis of the convergence for the 
eigenvalue problems by the AFEM. 

The purpose of this paper is to propose and analyze a type of AFEM to solve the 
eigenvalue problems based on the recent work on the multilevel correction method 
(see [3D]) and the two-grid correction method (see [IB])- In this new scheme, the 
cost of solving eigenvalue problems is almost the same as solving the associated 
boundary value problems. Here, we adopt the techniques in [201 EH EE] to prove the 
convergence and optimal complexity of the new AFEM for the eigenvalue problems. 
Our analysis is also based on the relationship between the finite element eigenvalue 
approximation and the associated boundary value problem approximation (c.f. [2111 

EZ])- 

The rest of the paper is arranged as follows. In Section 2, we shall describe some 
basic notation and the AFEM for the second order elliptic problems. In Section 3, 
we introduce a type of AFEM for the second order elliptic eigenvalue problems based 
on the multilevel correction scheme. The convergence analysis of this type of AFEM 
for eigenvalue problems will be given in Section 4 and Section 5 is devoted to proving 
the corresponding optimal complexity. In Section 6, some numerical experiments 
are presented to test the theoretical analysis. Finally, some concluding remarks are 
given in the last section. 

2 Preliminaries 

In this section, we introduce some basic notation and some useful results of AFEM 
for the second order elliptic boundary value problem. 

Let Q C 1Z d (d > 1) denotes a polytopic bounded domain. We use the standard 
notation for Sobolev space W S,P (Q) and their associated norms and seminorms (see 
e.g, DQ EH]). We denote H S ({1) = W s > 2 (n) and = {v G if 1 (£2) : v\ dn = 0}, 

where v\qq is understood in the sense of trace, |H| s ,n = IM|s,2,n and |H|o,n = 
II ^ 1 1 0,2,0- For simplicity, following Xu [46J, we use the symbol < in this paper. The 
notation A < B means that A < CB for some constant C independent of the mesh 
sizes. Throughout this paper, we shall use C to denote a generic positive constant, 
independent of the mesh sizes, which may varies at its different occurrences. We con- 
sider finite element discretization on the shape regular family of nested conforming 
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meshes {%,} over Q: there exists a constant 7* such that 

S 7 * Vi fc 1 

ft 



— < 7 * vtgMt,, 

PT W 



where /i^ denotes the diameter of T for each T 6 71, and is the diameter of the 
biggest ball contained in T, h := max{/i T : T G %,}■ In this paper, we use Eh to 
denote the set of interior faces (edges or sides) of Th- 

The following lemma is a result of Sobolev trace theorem (e.g. [U I3"T]). 

Lemma 2.1. If s > 1/2, then 

\\v\\o,dT < h T 1/2 \\v\\ , T + h s i 1/2 \\v\\ s , T \/v e H S (T), T G T h . (2.1) 
For / G L 2 (Q), we define the data oscillation (see, e.g., [331 [35]) by 



osc(/,Tft) := ||/»(/-A)lkn=(X)lM/-/T)llo,T) 

\T&r h J 



1/2 



(2.2) 



where fh denotes a piecewise polynomial approximation of / over Th and = fh\r- 
We will denote fx be the L 2 projection of / onto polynomials of some degree, which 
leads to the following inequality (see [20]): 

osc(/! + h, %) < osc(A, %) + osc(/ 2 , T h ) V/i, h e ^ 2 (fi). (2.3) 

2.1 Short survey on linear elliptic problem with adaptive 
method 

In this subsection, we shall present some basic results of AFEM for the second order 
elliptic boundary value problem. Here, for simplicity, we consider the homogeneous 
boundary value problem: 

Lu := — V(A ■ Vu) + <pu = f in Cl, , , 

u = ondfl, 1 ' 

where A = (aij)dxd is a symmetric positive definite matrix with G W 1,DO (n) = 
1,2- ■ -tZ), and < <p G 

The weak form of (12.41) is: Find w G Hq(Q) such that 

o(«, «) = (/,«) v.e^n), (2.5) 
where the bounded bilinear form is defined by 



a(u, v 



)= / (Vw- Vv + <puv)dQ. (2.6) 
Jn 



From the properties of A and ip, the bilinear form a(-, •) is bounded over Hq(Q) 
\a(w,v)\ < C a \\w\\ a! a\\v\\ a; n Vw,v G ff^fi), 

and satisfies 

c a |Mkn < \\w\\a,u < C a \\w\\ hn , 

where the energy norm || • || aj n is defined by ||w|| a ,n = a/ a(w, w), c a and C a are 
positive constants. From these properties, it is well known that ( 12.51) has a unique 
solution u G Hl(Vt) for any / G H~ 1 (Q). 

Let Vh C Hq(Q) be the corresponding family of nested finite element spaces of 
continuous piecewise polynomials over l~h of fixed degree m > 1, which vanish on 
the boundary of Q, and are equipped with the same norm || • || 0) n of space Hq(Q). 

Based on the finite element space Vh, we define the finite element scheme for (12.41) : 
Find Uh G Vh such that 

a(u h , v h ) = (f, v h ) \/v h eV h . (2.7) 

Define the Galerkin projection R h : Hq(Q) —> Vh by 

a{u-R h u,v h ) = Vv h eV h , (2.8) 

then we have Uh = RhU and 

\\Rhu\\ a ,u < \\u\\ a ,n Vu G H*(tt). (2.9) 

From (12. 9p . it is easy to get the global a priori error estimates for the finite element 
approximation based on the approximate properties of the finite element space Vh 

(c.f. pang). 

In order to simplify the notation, we introduce the quantity r) a (h) as follows: 

7] a (h) = sup inf WL^f - v h \\ a ,n. 

feLHn),\\fhn=i VheVh 

From [31 [HI [37], it is known that T] a (h) — > as h — > and the following propositions 
hold. 

Proposition 2.1. 

||(J - i4)L-7lUn < Va(h)\\f\\ 0> n V/ G L 2 (fi) (2.10) 

and 

\\u - R h u\\ 0:n < Va (h)\\u - R h u\\ atQ Vn G H*(n). (2.11) 
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Next we follow the classic routine to define the a posteriori error estimator for 
finite element problem (12. 7p . Let us define the element residual TZxiuh) and the 
jump residual J E {uh) by 

n T {u h ):=f-Lu h = f + V-{AVu h )-^u h in T e T h , (2.12) 
Je{u h ) := -AVul ■ u+ - AVul ■ v~ := [[AVu h ]] E -v E on Ee £ h , (2.13) 

where E is the common side of elements T + and T~ with outward normals v + and 
v~ , v E = v~ , and oj e '■= T + PI T~ that share the same side E. 

For the element T G 7/,, we define the local error indicator fj h (uh,T) by 

Vh(u h ,T) := (h 2 T \\n T (u h )\\l T + Yl hMlMuhWo,*) , (2-14) 
and the error indicator for a subdomain u C f2 by 

J 2 

rj h [u h ,u) : 



✓ \ 1/2 

E ^K' T ) ■ ( 2 - 15 ) 

\TeT h ,Tcw J 



Thus rjh(uh, fi) denotes the error estimator of with respect to 7^. 

Now we summarize the reliability and efficiency of the a posterior error estimator 
(see, e.g., [331 EH |H]): 

Lemma 2.2. ([20]) When ho is small enough, there exist mesh independent con- 
stants Ci, 6*2, C3 such that 

\\u-u h \\ at n < C^hiu^tt), (2.16) 

and , for any T e % 



C 2 2 f,l{u h ,T) - Cl Y, h 2 T \\TZ T (u h )-n T (u h )\\l T < \\u-u h \\l„ T , (2.17) 

TET h ,TCu T 

where ojt contains all the elements that share at least a side with T, C\,C2,Cz > 
depends only on the shape regularity 7*, C a and c a , lo is the L 2 -projection of w onto 
polynomials of degree m on T. 

As a consequence of \2.17\) , we have 



Cif! 2 h (u hl fi) - C 3 2 osc(/ - L(u h ), Th) 2 < \\u - u h \\ 2 (2.18) 



where C 2 = tt and Cq is a constant depending on the shape regularity of the mesh 

- Co 

lh- 
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It is obvious that if the right hand side term / of ( 12. 4 j) is a piecewise polynomial 
function over Th, (I2.18P can be simplified to 

cifj 2 h (u h ,n) -cio SC (L(u h ),r h ) 2 < \\u-u h \\l n . (2.19) 

There are adaptive algorithms in [TU [2TJ, [331 ES] to solve (12. 7p . which introduce 
two type of marking strategies to promise reduction of both error and oscillation. 
For these two type of methods, both convergence and optimal complexity of the 
adaptive finite element algorithm have been obtained (see, e.g., [2~T | [33| l35 | l4TJ | l4~T] ). 
However, the oscillation marking is not necessary which has been proved by Cascon 
et al. |l4j . Thus, the adaptive algorithm without oscillation marking which is 
adopted in this paper can be stated as follows (c.f. [TU [2D]). 



Adaptive Algorithm C 

Choose parameter < 9 < 1: 

1. Let k = 0, pick an initial mesh Th and start the loop. 

2. On the mesh %, k , solve the problem (12. 7p for the discrete solution u\ lk . 

3. Compute the local indicators fjh k (uh,T). 

4. Construct the submesh Th k C Th k by Marking Strategy Eq with parameters 
9. 

5. Refine 7~h k to generate a new conforming mesh 7~h k+1 by procedure REFINE. 

6. Let k = k + 1 and go to step 2. 



As in [13], the procedure REFINE used in Adaptive Algorithm Co is not 

required to satisfy the Interior Node Property of [331 HE]. Here we use the iterative 
or recursive bisection (see, e.g., [321 H2]) of elements with the minimal refinement 
condition in the procedure REFINE. The Marking Strategy adopted in Adap- 
tive Algorithm Co was introduced by Dorfler [2JJ and Morin et al. [33] and can 
be defined as follows. 



Marking Strategy E 

Given parameter < 9 < 1: 

1. Construct a minimal subset Th from Th by selecting some elements in Th 
such that 

^2 VH( u h,T) > 9f) 2 H (u h ,tt). 

T&f H 

2. Mark all the elements in Th- 
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In order to analyze the convergence of the AFEM, we need the following lemma. 

Lemma 2.3. ([20]) There exits a constant C only depending on the equation pa- 
rameters and the mesh regularity 7* such that 

osc(L(v h ),r h )<C L \\v h \\ ajQ Vv h eV h . (2.20) 

The convergence of Adaptive Algorithm Co has been proved by Cascon et al 
and can be stated as follows. 



Theorem 2.1. ([Tffl ) Let {uh k }k&N be a sequence finite element solutions of fli?.^| ) 
based on the sequence of nested meshes {Th k }k<m Q produced by Adaptive Algorithm 
Cq. Then, there exist constants 7 and £ G (0, 1), depending on the shape regularity 
of meshes, the data and the parameters used in Adaptive Algorithm C , such 
that any two consecutive iterates k and k + 1 have the property 



\u 



u h k+1 \\ln + lvl k+1 (uh k+ ^ Q ) ^ t 2 (\\ u - u h k \\ln + ivl k (u hk ,n)^, (2.21) 



where Nq = {0, 1, 2, 3, • ■ ■ } and the constant 7 has the following form 



with some constant S G (0, 1). 



3 The eigenvalue problem and adaptive finite el- 
ement method based on multilevel correction 

In this section, we introduce a type of AFEM based on multilevel correction scheme 
for the linear second order elliptic eigenvalue problem. 

We are concerned with the following eigenvalue problem 

— V(A • Vm) + ipu = Xu in fl, 

u = on dn, (3.1) 
f n u*dn = 1. 

The corresponding weak form can be written as: Find (X,u) G TZ x Hq(Q) such 
that ||u||o,fi = 1 and 

a(u,v) = \{u,v) Vi> G Hq(Q). (3.2) 
As we know the eigenvalue problem (13. 2p has a countable sequence of real eigenvalues 

< Ai < A 2 < A 3 < ■ ■ • 
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and corresponding orthogonal eigenfunctions 

Ui,u 2 ,u 3 , ■■■ , 

which satisfy (ui, Uj) = o~ij,i,j = 1, 2, • • • . 

Now we state an useful Rayleigh quotient expansion of the eigenvalue which is 
expressed by the eigenfunction approximation (see [31 1301 I4"8]). 

Lemma 3.1. Let (A, u) be an eigenpair of 113. || ). Then for any w G (Q), we have 

a(w,w) ^ a(w — u,w — u) X(w — u,w — u) (3 3) 

(w,w) (w,w) (w,w) 

The standard finite element discretization for ( 13. 2 p is: Find (\h,Uh) S TZ x Vh 
such that \\uh\\o = 1 and 

a(u h ,v h ) = X h (u h ,v h ) Wv h eV h . (3.4) 

We can also order the eigenvalues of (13. 4p as an increasing sequence 

< Xi th < X 2 ,h < X 3ih < ■ ■ ■ < X nhjh , n h = dim!4, 

and the corresponding orthogonal eigenfunctions 

Ul t h, U 2 ,h, u 3,h, ■ ■ ■ , U nh) h 

satisfying (u ijh , u jyh ) = i, j = 1, 2, • • • , n h . 

From the minimum- maximum principle (see [31 [T5]) and Lemma 13.14 we have 

Xi < X i)h < Xi + Ci\\ui - u iih \\l 7n , i = 1,2, , • • • ,n h (3.5) 

with C{ constants independent of mesh size h. 

3.1 Adaptive multilevel correction algorithm for eigenvalue 
problem 

The adaptive procedure consists of loops of the form 

Solve — > Estimate — > Mark — > Refine. 

Similarly to Marking Strategy E , we define Marking Strategy E for (13 .4p 

to enforce the error reduction as follows: 
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Marking Strategy E 

Given a parameter < 9 < 1: 

1. Construct a minimal subset Th from Th by selecting some elements in Th 
such that 

VH(uh,T) > 9r} 2 H (u h ,Q), 

Te? H 

where r]h{uh,T) and r]h(uh,Q) denote the error indicator of the eigenfunction 
approximation on T and Q, respectively. 

2. Mark all the elements in Th- 



Then we present a type of AFEM to compute the eigenvalue problem in the 
multilevel correction framework. 



Adaptive Algorithm C 

1. Pick up an initial mesh Th with mesh size h . 

2. Construct the finite element space Vh and solve the following eigenvalue 
problem to get the discrete solution (X ho ,Uh ) ElZx Vh such that ||tth ||o,n — 1 
and 

a(u ho ,v ho ) = \ ho (u ho ,v ho ) Vv ho eV ho . (3.6) 

3. Let k = 0. 

4. Compute the local error indicators r] hk (u hk ,T). 

5. Construct Th k C 7h fe by Marking Strategy E and parameter 9. 

6. Refine 7~h k to get a new conforming mesh l~h k+1 by procedure Refine. 

7. Solve the following source problem on %, k+1 for the discrete solution Uh k+1 € 

a (wft fc+ i. u h fc+ i) = VKfc.Vftfc+i) (3.7) 

8. Construct the new finite element space Vh ,h k+1 = Vh +span{M/ lfc+1 } and solve 
the eigenvalue problem to get the solution (\ hk+1 ,Uh k+1 ) ElZx Vh (h h k+1 such that 
IK fc+ illo,n = 1 and 

a ( u h k+1 , v h h0:hk+1 ) = ^h k+1 (u hk+1 ,Vh ,h k+1 ) Vv h0jhk+1 G V h(hhk+1 . (3.8) 

9. Let k = k + 1 and go to Step 4. 
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Local error indicator i]h(uh k ,T) in Adaptive Algorithm C will be given in the 
next subsection. For the aim of error estimate, we define 

M(Ai) = {weH^Q) : w is an eigenfunction of (13. ip 

corresponding to the eigenvalue Aj} 



and the quantity 



8 h (Ai) = sup inf - VftHo.n- 

«>eM(A l ),|k||o,n=l ,; ' ieVh 



In the following analysis, we only need some crude priori error estimates stated 
as follows. 

Lemma 3.2. The obtained eigenpair approximation (A^ fc , Uh k ) (k = 0, 1, ■ ■ ■ ) after 
each adaptive step in Adaptive Algorithm C has the error estimate 

\\u-u hk \\a,n < Sho(X), (3.9) 
||«-«/nllo,n ^ Va(h )\\u - u hk \\ a>Q , (3.10) 
l^-VJ ~ llw-w/ijlln- ( 3 - n ) 

Proof. Based on the error estimate theory of eigenvalue problem by finite element 
method (c.f. [31 H]), the eigenfunction approximation of problem (13. 6p or (13. 8p has 
the following estimates 

\\u-u hk \\ at n < sup inf \\w - v h \\ aj n 

< sup inf \\w -v h \\ at n <5 ho (\). (3.12) 

weM(X) v h£V hQ 

and 

\\u-u hk \\ Qin < rj a (h )\\u - u hk \\ a ,n, (3.13) 

where 

Vaih) = sup inf \\L~ l f - v\\ a ,n < Va(ho). (3.14) 

fev,\\f\\ , n =i v ^ v h ,h k 

From (I3TT2|) . (I3TT3|) . and (ETHjl . we can obtain fl33]) and (IBTTOjl . The estimate ( BHD 
can be derived by Lemma 13.11 and (I3.10p . □ 

3.2 A posteriori error estimate for eigenvalue problem 

Now, we are going to give an a posteriori error estimator for the eigenvalue prob- 
lem. The a posteriori error estimators have been studied extensively (see, e.g., 
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mi |23l EH ESI |2Sl EI]). Here we use the similar way in [2Q] to derive the a 
posteriori error estimator for the eigenvalue problem by Adaptive Algorithm C 
from a relationship between the elliptic eigenvalue approximation and the associated 
boundary value approximation. 

In this paper, we set H = and h = hk+i- Let K : L 2 (Q) — > Hq(Q) be the 
operator defined by 

a{Kw, v) = (w, v) Vw,v e H^tt). (3.15) 

Then the eigenvalue problems (13. 2p and (13.41) can be written as 

u = XKu, u h = \ h R h Ku h . (3.16) 

In the step 7 of Adaptive Algorithm C, we can view (13. 7p as the finite element 
approximation of the problem 

a(w H ,v) = (X H u H ,v) VuGi^O). (3.17) 

Thus we have w H = XhKuh and 

u h = R h w H . (3.18) 
Similarly we can also define w h as 

a(w\v) = (\ h u h ,v) VveH^n). (3.19) 

Let r(ho) '■— f] a (h ) + \\u — uh\\o,ci- Obviously from Lemma I3T21 we know r(h ) 1 
if ho is small enough. 

Theorem 3.1. We have the following estimate 

\\u - u h \\ at Q = \\w h - R h w h \\ at n + O(r(h ))(\\u - u H \\ a! n + \\u - Ufc||a,n)- (3.20) 

Proof. From the definition (13.81) . we have 

it — Uh = u-w h + w h - R h w h + R h w h - R h w H 
+R h w H - u h 
= u-w h + w h -R h w h + R h (w h -w H ) 

+u h - u h - (3-21) 

The following equality holds 

u - w h = XKu - \ h Ku h = \K(u - u h ) + (A - X h )Ku h , (3.22) 

which together with the fact \\K(u — Uh)\\ a ,n ^ — i*^||o,q, (13. lOj) . and (I3.1ip leads 
to 

\\u-w h \\ ai n < Cr(h )\\u - u h \\ a ,a- (3.23) 
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Similarly, we have 

\\R h w h - RhW H \\a,n < \\w h - w H \\ at n = \\\hKu h - \ H Ku H \\a,n 

< Cr(h )(\\u-u h \\ a>n + \\u-u H \\ a> n). (3-24) 

Since Uh ~ Uh E Vh ,h, fl3.10j) . and (13. lip , the following inequality holds 

a{u h - u h , u h - u h ) = (XhUh - XhU H , u h - u h ) 

< \K - ^H\(u h ,U h - U h ) + \ H {u h - U H ,U h - U h ) 

< Cr(h Q )(\\u - u h \\a,a + \\u - u H \\ a ,n)\\u h - u h \\ a ,a- (3.25) 
Thus with the coercivity of a(-,-), we have 

IK - Wfc||o,n < Cr(h )(\\u - u h \\ a ,n + ||« - whIUo)- (3-26) 

Finally from <K21h . (I3T231 . ([321, and <K2^i . the desired result fEOUD can be ob- 
tained and the proof is complete. □ 

Theorem 13.11 builds a basic relationship between \\u — Uh\\ and \\w h — Rhiu h \\, the 
former is the error between the ture and the discrete eigenf unctions, while the latter 
the error between w h and its finite element projection, which has been well analyzed. 
Since the difference between \\u — Uh\\ and \\w — RhW h \\ is a higher order term, as 
in [20], we follow the procedure of the analysis of convergence and complexity for 
the source problem. 

We define the element residual TZr{uh) and the jump residual ^(w/i) as follow: 

K T {u h ) := X h u h - L(u h ) = XhUh + V ■ {AVu h ) - ipu h in T E %, (3.27) 
Je{u h ) := -AVut ■ v + - AVul ■ v- := [[AVu h }} E ■ v E on E E £ h , (3.28) 

where E, z/ + and v~ are defined as those of Sect. 12.11 

For each element T E Th, we define the local error indicator r]h(uh,T) by 

V 2 h (u h ,T) := h 2 T \\TZ T (u h )\\l T + h E\\J E (u h )\\l E . (3.29) 

E£S h ,EcdT 

Then on a subset u C f2, we define the error estimator r] h (uh,oo) by 



r) h (u h ,uj) := 




(3.30) 



As same as (12.31) and (I2.20p . we have the similar inequalities of the indicator i]h(vh, oj) 
for any v h EV h . 

Lemma 3.3. The following inequalities for the indicator rjh{vh,uj) hold 

r] h (w h + v h , u) < r} h (w h , u) + r}h{vh, u) Vw h E V h , v h E V h , (3.31) 
r) h (v h ,Q) < C R \\v h \\ at n \/v h EV h . (3.32) 
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Proof. The first inequality (I3.3ip can be obtained from the definition of 77^. Now we 
prove the second inequality ( I3.32p . 

It is obvious that the inverse estimate implies 

\\L(v h )\\ ,T < C A h^\\Vv h \\ ,T + C c \\v h \\ , T VTeT h , (3.33) 

where Ca depends on A and the shape regularity constant 7*, C c depends on the 
coefficient <p. Namely, there exist some constants Ct and Cr depending on Ca and 
C c such that 

£ h T\\L(v h )\\l T < E & t\K\\1,t < C 2 R H\\lu- (3-34) 

T£T h TeT h 

From the trace inequality (12. ip and the inverse estimate, we have 

h E \\[[AVv h ]} E ■ v E \\\ E < C E \\AVv h \\l WE +C E h 2 T+ \\AVv h \\l T+ 

+C E h 2 T -\\AVv h \\l T - 
< ClII^H^, (3.35) 

where u E := T + U T~ denotes the patch including the elements sharing the edge E 
and the constant C E depends on A and the shape regularity constant 7*. 

Hence the desired result (13.321) can be obtained from (I3.29p . (I3.34p . and (I3.35P 
and the proof is complete. □ 

Based on Theorem 13. II and Lemma [3.31 the error estimate rj^Uh, SI) has the fol- 
lowing properties. 

Theorem 3.2. Let ho be small enough and h £ (0, ho]. Then there are mesh inde- 
pendent constants such that 

\\u - u h \\ a ,Q < Ctfh^ty + O(r(h ))\\u - u H \\ at n, (3.36) 

and 

C 2 V 2 (u h , n) - C 3 2 osc(LK), T h ) 2 <\\u- u h \\l n + O(r 2 (h ))\\u - u H \\l n . (3.37) 
Consequently we have 

|A - X h \ < vl(u h , n) + O(r 2 (h ))\\u - u H \\l Q , (3.38) 

and 

V 2 (u h , SI) - osc(LK), T h ) 2 < I A - A fc | + O(r 2 (h )) \\u - u H \\l Q . (3.39) 
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Proof. From flgjg) , (E2D}, (I53T]) . ( M) , (BHD , and (13T32]) . we have 

\\w h - R h w h \\ a , n < C 1 r] h (R h w h , Q) = Ctfhfah + R h w h - u h , fi) 

< Cxrjhiuhjfl) + Cir] h (R h w h - u h ,Q) 

< CiT} h (u h ,n) + C 1 C R \\R h w h -u h \\ a>n 

< CtfhfahtSi) +C 1 C R r(h )(\\u-u h \\ ai n + \\u -u H \\ a ,n), (3-40) 
which together with (I3.20p leads to the desired result (I3.36P with 

d = C 1 . (3.41) 
Combing m , qOS|) , (JCT]), and $£32} leads to 

(7 2 2 ^(^,^)-C 3 2 osc(LK),r ft ) 2 

< C 2 fj 2 h (R h w h , n) - C 2 osc(L(R h w h ), T h ) 2 
HC 2 L + C 2 R )(C 2 + C 2 3 )\\u h -R h w% n 

< \\w h -R h w h \\ln + C\\u h -R h w h \\l n 

< 2\\u - u h \\l n + Cr 2 (h )(\\u - u h \\l n + \\u - u H \\l n ), (3.42) 

where the constant C depends on C, Cl, Cr, C 2 and C 3 . From ( I3.42p . the desired 
result f )3.37p can obtained with 

C 2 C 2 
C2 = 2 + Cr*(h y Cl = 2~Tcks)' ^ 

From Lemma 13.11 and (I3.36p . we obtain (I3.38P and f!3.39j) can be derived from 
(E3D and (J33H). □ 

4 Convergence of adaptive finite element method 
for eigenvalue problem 

In this section, we give the convergence analysis of the Adaptive Algorithm C 
for the eigenvalue problem. 

Before establishing the error reduction of the Adaptive Algorithm C for the 
eigenvalue problem, we give some preparations. 



Similarly to Theorem I3.1[ we also give some relationships between two level ap- 
proximations, which will be used in the following analysis. 

Lemma 4.1. Let h,H G (0, ho]. If w h = XhKuh, w H = XhKuh, we have 
||w-Mfc||o,n = II w H - R h w H II O o 
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+O(r{h ))(\\u-u h \\ a> n + \\u-u H \\ a>a ), (4.1) 
osc(L(u h ),Th) = osc(L(R h w H ),T h ) 

+O(r(h ))(\\u - u H \\ a:Q + ||u - Ufc||o,n), (4.2) 

and 

Vh{u h ,n) = fj h {R h w H ,Q) + O(r(h ))(\\u-u h \\ at n + \\u - u H \\ a>n ). (4.3) 

Proof. First we have 

u - u h = u - w H + w H - R h w H + R h w H - u h . (4.4) 

Similarly to (j3.23p . the following inequality holds 

||« - w H \\ a ^ = \\u - K\ H u H \\an = \\KXu - K\ H u H \\ a: n 
< \\Xu - X h uh\\o,q = O(r(h ))\\u - u H \\ a ,u (4.5) 

Combining fl3T26|) . fOj) . and gj| leads to (Q|l . 
By the property (12.201) of oscillation, we have 

osc(L(R h w H - u h ),Th) < C L \\R h w H - u h \\ a> n, 

which together with (I3.26P and the fact = RhW H implies 

osc(L(R h w H - u h ),Th) < r(h )(\\u - u H \\ a ,n + \\u - u h \\ a> n). (4.6) 

Hence from (j2T2Dl . (Oil , and 

L(u h ) = L(R h w h ) + L(R h w H -u h ), (4.7) 

we can obtain the desired result (14. 2p . 

Now we come to consider the relation (14.31) . Using (13.261) . (13.321) . and the fact 
Uh = RhW H , we obtain 

fj h (R h w H -u h ,tt) < r(ho){\\u-u H \\a,a + \\u-Uh\\ a fl)- ( 4 -8) 
Combining (I33T|) . (Q)l and the fact 

t7h(iih, J]) = r] h (R h w H + u h - R h w H , Q) < rj h (R h w H , SI) + r] h (u h - R h w H , SI) 
< rj h (R h w H , Q) + O(r(h ))(\\u - u h \\ a ,u + \\u - u H \\a,u) 

leads to the desired result (14. 3 p and the proof is complete. □ 

Now we are at the position to give the error reduction of Adaptive Algorithm 
C for the eigenvalue computations. 
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Theorem 4.1. For the successive eigenfunction approximations uh and Uh produced 
by Adaptive Algorithm C, there exist constants 7 > 0, «o, and a G (0,1), 
depending only on the shape regularity of meshes, C a , c a and the parameter 6 used 
by Adaptive Algorithm C , such that 

\\v>-u h \\l n + jr)l(u h ,n) < a 2 (\\u-u H \\ 2 aQ + ^r] 2 H (u H ,Q)) 

+a 2 l r 2 (h )\\u - u jH _J 0>n , (4.9) 

provided h ^ 1. 

Proof. Since w h = X^Ku^ and w H = XhKuh, we conclude from Theorem 12.11 that 
there exists constant 7 > and £ G (0, 1) such that 



IK-iWlin+ 7^4^,0) < e(\\w H -R H w«\\l n 

+^ffj 2 H (R H w H ,Q)). (4.10) 

From ( 14. ip and (14. 3p . there exists a constant C > such that 

11^-^11^ + 7^(^,0) 

< (1 + 5 1 )\\w H - iWllin + 6(1 + ^V^oXh - «*|lin + h~ "ffllln) 
+(1 + fc^i?^* O) + (77(1 + 5r V(M(II^ - ttfcHin + ||u - u H \\l^) 

< (1 + 50 - iWHlln + ivl(R h w H , O)) 

+C 4 5r 1 r 2 (/ i0 )(||« - Whll^n + ||« - ujfllin), (4.11) 

where C4 depends on the constants C and 7 and the Young inequality is used with 
S\ G (0, 1) satisfying 

(i+s 1 )e<i- 

The similar argument leads to 

\\w H - R H w H \\l n + if&{R H w H , n) 

< (1 + 5 2 )\\u - u H \\l n + (7(1 + 6^)r 2 (h )(\\u - u H \\l Q + \\u - u u _, |Q 
+(1 + 5 2 )^r, 2 H (u H , SI) + (77(1 + ^V^oXHw - u ff ||^ n + ||u - ujr.Jin) 

= (i + $ 2 )(|| u - u H \\ln + ir&(uH,n)) 

+C i 5 2 - 1 r 2 (h )(\\u-u H \\l n + \\u-u H _ 1 \\ 2 a ,n)i (4-12) 

where Uff^ denotes the eigenfunction approximation obtained on the mesh level 
7fo fc _ 1 before Th and <5 2 G (0, 1) satisfies 

(1 + 5 1 )(1 + S 2 + C^ l r 2 (h ))e < 1. (4.13) 
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Combing fl4TTUj) and (HTTj) leads to 

lk-^lla,n + 7^K>^) 

< (i + *i)e a (||«^ - R H w H \\ 2 a>n + if H (R H w H , ny 

+C^ 1 r 2 (h ){\\u - u h \\l n +\\u- u H \\l n ). (4.14) 
From ( I4.12P and ( I4.14j) . we have 

h-UhWln + lvKuh,^) < (l + 6i)?(p. + 62)(\\u-UH\\l a + mUuH,n)) 

+C 4 5 2 1 r 2 (h )(\\u - u H \\l } u + h - Wff-JaA)) 
+^r 1 r 2 (/ i0 ) (||t* - u h \\l n + \\u- u H \\l^ . (4.15) 

Consequently, 

(i - c^-V^o)) Ik - «k|lin + 7^K, n) 
< ((1 + 50(1 + 5 2 + C^-V^/iq))^ + C^V 2 ^)) \\u - u H \\l n 
+(l + 5 1 )(l + 5 2 )e 2 7^K,fi) 

+C 4 (1 + 6 1 )6^ 2 r 2 (h )\\u - u H _ x f a ^ (4.16) 



that is 



\\ u ~ u>h\\l,n + n jzi 2 ,, , vl(uh, ft) 

. (1 + 50(1 + 6 2 + Q^-V^q))^ + C4*rV 2 (M H |,2 
1 - C A b 1 r 2 (h ) 



1 - C 4 6^r 2 (h 



r (Mlu-Uir-Jtn- ( 4 - 17 ) 



Since /io «C 1 implies r(h ) <C 1, we have that the constant a defined by 
= ( (1 + ft)(l + + C^r 2 ^))^ + C^Mfep) \ 1/2 

satisfying a 6 (0, 1). Therefore 

\\ u ~ u h \\l,Q + ~ n jL 2(u , vl(uh, ft) < a 2 ( ||u - u H \\ln 

, (l + ^)(l + ^ 2 7 2 

^(l + ^a + ^ + ^^^^ + ^rV 2 ^) H ' J 
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C4(1 + WW|| ||2 r . 1Q x 

ll«-«ff-ilLn- ( 4 - 19 ) 



i - c^rVOo) 

If we choose 



7 := - r Zi 2(hV ( 4 - 2 °) 
1 - CVi H(fto) 



we arrive at (14. 9 p by using the fact 

(l + <Ji)(l + <5 2 )£ 2 7 



<7 



(1 + 50(1 + 5 2 + C^V^o)^ + C&^iho) 
and setting 

2 _ C,(l + 5 1 )5 2 1 e 

Hence the proof is complete. □ 

Based on Theorem 14.11 we can give the following error estimate for Adaptive 
Algorithm C. 

Theorem 4.2. Suppose (X,u) G 1ZxHq(Q) be a simple eigenpair of (Xh k ,Uh k ) 
be a sequence of finite element solutions produced by Adaptive Algorithm C. 

When ho is small enough, there exist constants (3 > and a G (0, 1), depending on 
the shape regularity of meshes and the parameter 9 , such that for any two consecutive 
iterates k and k + 1 

4 k+1 < rfd^ + ayiho^u-Uh^Wln, (4-21) 

and 

d 2 hk+i + (5 2 r\h )dl k < a 2 {d 2 hh + ^(h )dl k i ), (4.22) 

where d\ = \\u — Uh h \\a t n + IVh ( u hh-> Then, Adaptive Algorithm C con- 
verges with a linear rate a, i.e. the n-th iteration solution {Xh n ,Uh n ) of Adaptive 
Algorithm C has the following error estimates 

dl + ^ihM^ < C a 2n , (4.23) 
lA^-Al < e?\ (4.24) 

where C = ||u - u/J|o, n + 7^o( u Ao> n )- 

Proof. It is obvious that (I4.2ip can be derived directly from Theorem 14.11 Now we 
choose a and (3 such that 

a 2 -(3 2 r 2 (h ) = a 2 , 
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-2 o2 2 

a p = a . 

This equation leads to 

_ 2 a 2 + y/a A + Aalr 2 (h ) 2 2 «o 

a = and p = . 

2 a 2 + ^/a 4 + 4:a 2 r 2 (h ) 

As we know a < 1 provided that a < 1 and fro is small enough. Then (I4.22p can be 
obtained with the chosen constants a and j3. 

The estimates (I4.23P and (14.241) are natural results from (I4.22p and (13.31) and the 
proof is complete. □ 



5 Complexity analysis 

Due to Theorems 12.11 and 14.11 we are able to analyze the complexity of Adaptive 
Algorithm C for eigenvalue problem via the complexity result of the associated 
boundary value problems. 

In this section, we assume the initial mesh size ho is small enough such that 

r(h )\\u - Uh^Wln < \\u - u hk ||* n . (5.1) 

Then from Theorem 14 .![ we have the following error reduction property of Adaptive 
Algorithm C 

\W - M ^ + illa,^ + lVh k+1 ( u h k+1 , n) < a 2 (Ik - u hk \\l n + jT] 2 hk (u hk ,n)) (5.2) 

with a 2 = a 2 + a^riho). 

Based on this contraction result, we also give the complexity analysis with the 
similar way of [2] and [2D]. Let 7~h k {k > 0) be the sequence of conforming nested 
partitions generated by REFINE starting from Th with /i Q <l. We denote Th k „ 
a refinement of 7h k ( m general nonconforming), M.(Th k ) the set of elements of Th k 
that were refined in Th k - Let Ih k+1 '■ C(O) D Hq(Q) — > Vh k+1 satisfy 

I hk+1 v = v on T^M(T hk ) WveV hk+1 

and set 

v f V hk U ((/ - I hk+1 )V hk+1 ), if Th kt . is nonconforming, 
hk '* \ V hk+1 , if Th kt , is conforming. 

In our analysis, we also need the following result (see, e.g., [2D1 EEH GEE EH HI])- 

Lemma 5.1. (Complexity of REFINER LetTh k (k > 0) be a sequence of conform- 
ing nested partitions generated by REFINE starting from Th , M.(Th k J the set of 
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elements of l~h k which is marked for refinement and Th k , be the partition created 
by refinement of elements only in Ai(Th kt ). There exists a constant Co depending 
solely on Th such that 



#T hk+1 - Wk < C J2i#T h ^ - #T hi ). (5.3) 

i=0 

Here and hereafter in this paper, we use #T to denote the number of elements in 
the mesh T. 

In order to analyze the complexity of Adaptive Algorithm C, we first review 
some results related to the analysis of complexity for the boundary value problem 
(12. 5p . For the proofs, please read the papers [H] and 



Lemma 5.2. (fT$ ) Let Rh k u G Vu k and Rh k ,u G Vh k , be the discrete solutions 
of h2.5\) on the meshes l~h k and its refinement Th k , with marked element M.{J~h k J- 
Then we have 

\\R hk u-R hk u\\ln<Ci vl k (R hk u,T). (5.4) 

TeM(T hk J 

Lemma 5.3. (fl^ \2Dj) Under the same assumptions as in Lemma \5.2\ and the 

energy decrease property 



\\ u ~ Rh k ,M\a,a + 7oosc(/ - L(R hk t u)) 

< H(\\u - R hk u\\l n + 7 osc(/ - L(R hk u) f) (5.5) 

with 7o > and £q G (0, |). Then the set A4(J~h k J of marked elements satisfy the 
Dorfler property 

J2 vl k (Rk k u,T) > 6 J2 fh h (Rh k u,T), (5.6) 

T&M(T hk J TeT hk 

where ~ e = c o{ clHi+2Clcim mth Co = max{l, f }. 

As in the normal analysis of AFEM for boundary value problems, we introduce a 
function approximation class as follows 

At, := {v G #o(fi) : \v\ sn < oo}, 

where 7 > is a constant and 



f| Si7 = supe inf _ {^rf-^rf h 



>0 {rcr ho :inf(||«-« r ||?+(7+l)osc(i(« r ,r)) 2 ) 1 /2< £ } 



h f 



20 



and T C Th Q means T is a refinement of Th - From the definition, for 7 > 0, we see 
that At = A\ and we denote A s as A{, \v\ a as \v\ s ^ for simplicity. Hence the symbol 
A s is the class of functions that can be approximated within a given tolerance e by 
continuous piecewise polynomial functions over a partition T with the number of 
degrees of freedom #T — #7h < e~ l / s \v \l. 

In order to give the proof of optimal complexity of Adaptive Algorithm C 
for the eigenvalue problem ( 13. ip . we should give some preparations. Associated 
with the eigenpair approximation (Xh k ,Uh k ) of (13. 4p in the mesh Th k , we define 
w hk = K(X hk u hk ) as in (13Tm 

Using the assumption ( 15. ip and the similar procedure as in the proof of Theorem 
14.11 when ( I4.3p is replaced by (|4.2|) . we have 

Lemma 5.4. Let (Xh k ,Uh k ) G TZ x Vh k and (Xh k ,,Uh k J G TZ x V/^, 6e t/ie discrete 
solutions of 113. 1\) produced by Adaptive Algorithm C over a conforming mesh Th k 
and its (nonconforming) refinement Th k t with marked element Ai(Th k „). Supposing 
they satisfy the following property 

\\ u - u h k , t \\l,u + l*osc(L(u hk J, T hk J 2 

< & 2 (||« - u hk + 7 *osc(LKJ, T hk ) 2 ), (5.7) 

where 7* > ; /3* > are some constants. Then the associated boundary value 
problem approximations Rh k w hk and Rh k t w hk of w hk have the following contraction 
property 

\\w h " - R hk ^w h » ||* n + lif osc(L(R hk w hk ), % K J 2 

< g(||«; fc *-i^t£; fcfc ||2o + 7.08c(L(t*0,7^) a ) (5.8) 



with 



> , (1 + *i + C A 5^r\h ))Pl + C&Vjho) V /2 
' ^ 1 l-C 4 r(/io)(l + r(/io)) ; ' [ ' 



where the constant C4 depends on S± G (0, 1) and 7* as m £/ie proof of Theorem \4-l 
Proof. From (12.20p . ( I4.ip . and ( I4.2p . there exists a constant C > such that 

||^ - R hk y« {fa + ^osc(L(R hk ^w hk ), % k J 

< (1 + 5 x )\\u - u hkt + 6(1 + 6^)r 2 (h )(Wu - u hkt Win + II" - Uh h Win) 
+(1 + 5 1 ) 7 *osc(LK fe> J, VJ 2 + (7(1 + 5r 1 )r 2 (/i )(||« - ||^ 

+ lk-^JIa,n) 

< (1 + 5i) (||w - w hk> , || 2 >n + 7*osc(L(w^J, 7;j 2 ) 

+C^ 1 r 2 (h ) (\\u - u hk Wl,n + W u ~ «hjla,n) 
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< (1 + S 1 )f3 2 (\\u - UhJJn + 7*osc(L(w hfc ), T hk ) 2 ) 
+C^ 1 r\h )/3l(\\u-u hk \\l n + ^osc(L(u hk ) } T hk ) 2 ) 
+C 4 5^r 2 (h )\\u-u hk \\l n 

< ((1 + 5 1 + C^r^hoM + C 4 5^r 2 {h )) (\\u - u hk ||* n 

+ 7 ,osc(LKJ,r h J 2 ), (5.10) 



where C4 depends on the constants C and 7* and the Young inequality is used. 
Similarly from f[2T2Uj) . (Q]l . and ([SI]), we have 

ll« - ^ lla,n + 7*osc(L(u fc J, 7;j 2 

< \\w h * -R hh w h % u + % osc(L(R hk w h >°), T hk ) 2 

+C 4 r 2 (h ) (\\u - u hk \\ 2 a n + ||w - || 2 >n ) 

< \\w h * - R hk w h % n + ^o S c(L(R hk w h *),T hk ) 2 

+C 4 r(h )(l + r(h ))\\u-u hk \\l n . 

Then the following inequality holds 

h - u h k \\ln + l*osc(L(u hk ), T hk f 

+7*osc(L(#^),T ftfc ) 2 )- (5.11) 
Combining (I5.10p and (15. lip leads to 

||^ - i^VilU + 7*osc(L( J R, fe ^^)) 2 



(1 + *i + c^rV^oM + gA' V 2 (^o) 

l-<7 4 r(/i )(l + r(V)) 
+^osc(L(R hk w h "),T hk ) 2 ). (5.12) 

This is the desired result (15. 8p and the proof is complete. □ 

We present the following statement which is a direct consequence of Lemmas 15.31 
and I 



Corollary 5.1. Let (\ hk ,u hk ) 6 TZ x V hk and (\h k *,u hk J G TZ X V hk „ 6e as m 



Lemma 5J^. Suppose they satisfy the following decrease property 
\\ u - u h k J\l,n + 7*osc(L( M/lfe J,7; fci J 2 < f3 2 (\\u-u hk \\l n + %osc(L(u hk ),T hk ) 2 ), 

where the constants 7* > and @ 2 G (0,1/2). T/ien the set M.{Th k ,) of marked 
elements satisfies the following inequality 
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where the constant 9 = ^ t^tttt-, — T an d Co = max{l, — } with /3* and 7* which 
are the same as in A5.8\) and 115. 9\) with 5\ being chosen such that 01 G (0, 1/2). 

Lemma 5.5. (Upper Bound of DOF). Letu G A s andTh k be a conforming partition 
from Th - Let l~h k , be a mesh created from Th k by marking the set M.(Th k J according 

to Marking Strategy E with 9 G (0, ^(^+{1+2^ c a h) ) • Then we have 

-#% <C(||n-^ fc ||^ + 7 osc(LKJ,r,J 2 )" 1/(2s Vls /s , (5-14) 
where the constant C depends on the discrepancy between 9 and c i {c' 1 +{i-^2C i c i )^) 1 



Proof. We choose /3, 0i G (0, 1) such that 0i G (0,/3) and 

9 < —. ^ -(1 - (5 2 ). 

Ci(Cf+ (l + 2ClC?h 

Let 

Let 5i G (0, 1) and 5 2 G (0, 1) be the constants such that (14.131) and 

(1 + Sl f(l + Szfff < P\ (5.15) 

which implies 

(l + 5 1 )(l + 5 2 )/3 1 2 < 1. (5.16) 

Let Th e be a refinement of 7h with minimum degrees of freedom satisfying 

\W - u he H 2 ^ + (7 + l)osc(LKJ, T he f < e 2 , (5.17) 

where Uh e denotes the solution of eigenvalue problem f l3.4j) over the mesh 7/ le . By 
the definition of A s , we can get that 

m e -4T h0 < f-^ft) 7 (|| M -^J| 2 in + 7 osc(LK fc ),r,J 2 )" 1/(2s Vl s 1/s - (5.18) 

Let 7h fe . be the smallest (nonconforming) common refinement of Th k and Th E ■ 
Since both Th k and 7h e are refinements of 7h , the number of elements in 7~h k , that 
are not in 7/^ is less than the number of elements that must be added to go from 
T ho to T he , namely, 

#% k>+ - m k < - #V (5-19) 
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Let w he = K(\ hs UhJ- From definition, we can easily get 

osc(L(R hk>+ w h *),T hk J < osc(L(R he w h °),T hh<+ ) 

+osc(L(R hk+ w h ° - R he w h *),T hk J 
< osc{L{R he w h °) 7 T hkt+ ) 

+C L \\R hkt+ w h ° - R hs w h *\\ a , Q , (5.20) 



where Lemma 12.31 is used. Then by the Young inequality, we have 

osc(L(R hk+ w h ^,T hk y < 2osc(L(R he w h ^,T hk y 

+2Cl\\R hkt+ w h * - R he w h *\\l n . (5.21) 

Since Th k+ is a refinement of 7/i e , L 2 -projection error are monotone and the fol- 
lowing orthogonality 

\\w h * - R hki+ w h *\\l n = \\w h ° - R hB w h *\\l n - \\R hk , + w h * - i^^Hin (5-22) 

is valid, we arrive at 

\\w h * - R hh!+ w h % n + ^osc(L(R hK+ w h ^,T hk ^ 2 

L 

< \\w h ° - R he w h % Q + ^osc(L(R h£ w h ^),T h J 2 . (5.23) 

Note that (12.221) implies 7 < and we obtain that 

\\w hs ~ R hK+ w h % n + 70Sc(L(R hk+ w h '), T hh J 

< \\w h ° - R he w h *\\l n + -^osc(L(R hs w h c),T hE ) 2 

< \\w h ° - R he w h °\\l >n + (7 + a)osc{L{R he w h *), T h f (5.24) 

with er = — 7 G (0, 1). Applying the similar argument in the proof of Theorem 
14.11 when f 14 . 3 1) is replaced by (14. 2p . we then obtain 

\\ u ~ u h k , + \\l,n + losc{L{u hk>+ ), T hk + ) 

< Pl (\\u - u hE ||* n + (7 + <j)osc{L{u he ), T he ) 2 ^j 

< Po(h - u hs \\la + (7 + l)osc(LKJ,T,J 2 ), (5.25) 



where 



A) 



1/2 

(1 + ft) ((1 + 6 2) + C 5 5^r 2 (h )) + C 5 5^r 2 (h ) 
1-C 5 5^(h ) 
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and 

= 7 
7 l-C 5 5^r*(h o y 

with C5 the constant depending on Cl similar to C4 in the proof of Theorem 14. 11 
Combing (15.171) and (I5.25P leads to 

\\u - u hk + \\ 2 an + nf08c(L(u hhi+ ), Th k + ) <$(\\u- u hk ||2 n + 7 osc(L(^J, 7^J 2 ) 

with = -fepoPi. 

It is seen from h <C 1 and (I5.16P that f3 2 G (0, ~). Thus by Corollary 15. II we have 
that Th k . satisfies 

E ^cn 2 >^ E ^( T ) 2 ' 

where $ = c (4it+lclciny 7 = i-c^rHhoy ^ = max ^' f > and 



1/2 



From the definition of 7 (see (14.201) ). we obtain that 7 < 1. On the other hand 
we have C3 > 1 and hence Co = Consequently, we can write 6 as 6 = 

C%(l-2p) 
Cf(^+(1+2C£C?))' 

Since /i 1 and (I5.15p . we obtain that 7 > 7 and /3 G (0, 775/?) ■ Using ( 13.41 j) 
and ( I3.43p . we get 

- C* 2 (l-2/3 2 ) C\ „ o2 . 

>—T^2 ^(l-/?) 



(7f(f + (1 + 2C 2 C< 2 )) " C 2 (f + (1 + 2C 2 Ci 2 

(2 + Cr 2 (^ ))C 2 ^ 
(2 + Cr*(/*))Cf(? + (l + 2ClCl" 

> —r-i ^(1 -/3 2 ) 

C!(^ + (l + 2ClCf 

ch 



^(l-f3 2 )>9. (5.26) 

C 2 (C 2 + (1 + 2C 2 C 2 )7) 

Note that Marking Strategy E selects a minimum set M.(Th k ,) satisfying 

E ^ T ?> d E ^-cn 2 ' 
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which implies that the nonconforming partition Th k , satisfies 

#^fc,* ~~ ^ ~~ #7/1* — #7h e _ #7fto 

<(^/?i) 1/ ^lk-^J| 2 O) n+70sc(LKJ,r ftfe ) 2 )- 1/(2s VL 1/s - (5-27) 

This is the desired estimate ( I5.14p with an explicit dependence on the discrepancy 
between 9 and q^+ff^^^ via ft. □ 

We are now in the position to prove the optimal complexity of Adaptive Algo- 
rithm C which is stated in the following theorem. Please refer the papers [14J and 
[20] for the proof. 



Theorem 5.1. ffm Let (\,u) G U x (H^(Q) n A s ) be some simple eigen- 
pair of 113.1}) and {(Xh k ,Uh k )}ken be the sequence of finite element approximations 
corresponding to the sequence of pairs {1Z x Vh k }kmo produced by Adaptive Al- 
gorithm C satisfying 113.10}) . Then under the assumption 115. 1}) . the n-th iterate 
solution (\h„, u h n ) °f Adaptive Algorithm C satisfies the optimal bounds 

ll«-«*Jlin + 70flc(L(«O.7^) a < WT hn ~ #r, )" 2s , (5.28) 

|Afcn-A| < (#r, n -#T, )- 2s , (5.29) 

where the hidden constant depends on the exact eigenpair (A, u) and the discrepancy 



between 9 and ttttttttht^ttti 



C|(C 1 2 +(1+2C£C 1 2 ) 7 ) ' 



Proof. We give the proof in the same way in [20]. From (13.371) and (15. 2p . we have 
u h k \\l,n + ivl k ( u h k > n ) <C[\\u-u hk \\l n + -fosc{L{u hk ),r hk ) 2 ), (5.30) 



\u 



where C depends on C 2 , C 3 , 7, and r(h ). 
Combining ( I5.14p and (I5.30p leads to 

#K. - #% k < CC* (\\u- u hk \\\ n + jVh h (u hk , n)J \u\ l ' s . (5.31) 
From (15. 2p we have 

/ \-l/(2s) 

[\\u-u hk \\l Q + jr]l k {u hk ,n)J 

< a^' s (\\u - u hn \\l n + 1V 2 hn (u hn , fl)) " V(2S) • (5.32) 
Employing Lemma [5. 1\ ( I5.14p . and ( I5.32p . we can obtain 

71-1 

#7h„ - #7h < 2^ (#7^=,* _ 7&J 

k=0 
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" ' -I/O) 



- \ U \1 /S ^2 (\\ u ~ u h k \\ln + ivl k ( u h k i n ) / 

/ (u 112 2 t ^\~ 1/(2S \ ,1/«V> * 



fe=i 

Combining the fact a < 1 leads to 

/ \-l/(2s) 

#T hn -#T ho < [\\u-u hn \\l n + -fri 2 hn (u hn ,n)j \u\ 
Sice osc(L(u h ),T hn ) < r] hn (u hn ,n), we have 

#T hn -#T ho < (||w-«fe n ||a,n + 70sc(u /ln ,7^ 

This is the desired result (I5.28P and (15.291) can be derived from (15.281) and Lemma 
13.11 Then the proof is complete. □ 

6 Numerical experiments 

In this section, we present sone numerical examples of Adaptive Algorithm C for 
the second order elliptic eigenvalue problems by the linear finite element method. 
Example 1. In this example, we consider the following eigenvalue problem (see 

!) 

-|Am + ||x| 2 m = Am in Q, 

u = on dtt, (6.1) 
\\ u \\o,n = 1, 

where f2 = 1Z 2 and \x\ = \Jx\ + x\. The first eigenvalue of (16. ip is A = 1 and 
the associated eigenfunction is u = /te - ^' I 2 with any nonzero constant k. In our 
computation, we set Q = (—5,5) x (—5,5). 

First, we investigate the numerical results for the first eigenvalue approximations. 
We give the numerical results for the eigenpair approximation by Adaptive Algo- 
rithm C with the parameter 9 = 0.4. Figure Q] shows the initial triangulation and 
the triangulation after 14 adaptive iterations. Figure |5] gives the corresponding nu- 
merical results for the first 19 adaptive iterations. In order to show the efficiency of 
Adaptive Algorithm C more clearly, we compare the results with those obtained 
with direct AFEM. It is observed from Figures [2], the approximations of eigenvalue 
as well as eigenfunction approximations have the optimal convergence rate which 
coincides with our theory. 

Example 2. In the second example, we consider the Laplace eigenvalue problem 
on the L-shape domain 

—Au = Xu in Q, 

u = on dQ, (6.2) 
Mlo.n = 1, 
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Initial mesh Mesh after 14 iterations 




Figure 1: The initial triangulation and the one after 14 adaptive iterations for 
Example 1 



where Q = (—1, 1) x (—1, 1)\[0, 1) x (—1, 0]. Since Q has a reentrant corner, eigen- 
functions with singularities are expected. The convergence order for eigenvalue 
approximations is less than 2 by the linear finite element method which is the order 
predicted by the theory for regular eigenfunctions. 

First, we investigate the numerical results for the first eigenvalue approximations. 
Since the exact eigenvalue is not known, we choose an adequately accurate approx- 
imation A = 9.6397238440219 as the exact first eigenvalue for our numerical tests. 
We give the numerical results for the first eigenpair approximation of Adaptive 
Algorithm C with the parameter 9 = 0.4. Figure [3] shows the initial triangulation 
and the triangulation after 12 adaptive iterations. Figure H] gives the corresponding 
numerical results for the first 20 adaptive iterations. In order to show the effi- 
ciency of Adaptive Algorithm C more clearly, we compare the results with those 
obtained by direct AFEM. We also test Adaptive Algorithm C for 5 smallest 
eigenvalue approximations and their associated eigenfunction approximations. Fig- 
ure [5] shows the corresponding a posteriori error estimator rjh(uh,Q) produce by 
Adaptive Algorithm C and direct AFEM. 

From Figures H] and |5j we can find the approximations of eigenvalues as well as 
eigenfunctions have optimal convergence rate which coincides with our theory. 

Example 3. In this example, we consider the following second order elliptic 
eigenvalue problem 

—V • (AVu) + ipu = \u in Q, 

u = onffl, (6.3) 
\\ u \\o,n = 1, 
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Eigenvalue errors Eigenfunction errors 




10 1 10 2 10 3 10 4 10 5 10 1 10 2 10 3 10 4 10 5 

Number of elements Number of elements 

Figure 2: The errors of the eigenvalue and the associated eigenfunction approxima- 
tions by Adaptive Algorithm C and direct AFEM for Example 1 

with 

1 + (xi - |) 2 (xi - §)(x 2 -\)\ 
{x x -\){x 2 -\) l + (x 2 -i) 2 )> 

(p = e (*i-|)(*a-|) anc i q = i) x i)\[ 0) i) x (-1, 0]. 

We first investigate the numerical results for the first eigenvalue approximations. 
Since the exact eigenvalue is not known neither, we choose an adequately accurate 
approximation A = 13.58258211870407 as the exact eigenvalue for our numerical 
tests. We give the numerical results for the first eigenpair approximation by Adap- 
tive Algorithm C with the parameter 9 = 0.4. Figure |6] shows the initial tri- 
angulation and the triangulation after 12 adaptive iterations. Figure [7] gives the 
corresponding numerical results for the first 18 adaptive iterations. In order to show 
the efficiency of Adaptive Algorithm C more clearly, we compare the results with 
those obtained with direct AFEM. 

We also test Adaptive Algorithm C for 5 smallest eigenvalue approximations 
and their associated eigenfunction approximations. Figure [S] shows the a posteriori 
error estimator f]h{ u hi produced by Adaptive Algorithm C and direct AFEM. 
From Figures [7| and El we can find the approximations of eigenvalues as well as 
eigenfunctions have optimal convergence rate. 




7 Concluding remarks 

In this paper, we present a type of AFEM for eigenvalue problem based on multilevel 
correction scheme. The convergence and optimal complexity have also been proved 
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Mesh after 12 iterations 




Figure 3: The initial triangulation and the one after 12 adaptive iterations for 
Example 2 

based on a relationship between the eigenvalue problem and the associated boundary 
value problem (see Theorem I3.ip . We also provide some numerical experiments to 
demonstrate the efficiency of the AFEM for eigenvalue problems. 
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